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g ! Abstract 

j^ _ In the vicinity of boundaries the bulk universaHty class of critical phenomena splits into several 

-)— > . 

"^ ■ boundary universality classes, depending upon whether the tendency to order in the boundary is 

-)— > 

rH . smaller or larger than in the bulk. For Ising universality class there are five different boundary 

F^ \ universality classes: periodic, antiperiodic, free, fixed and mixed (mixture of the last two). In this 

Q ' paper we present the new set of the universal amplitude ratios for the mixed boundary universality 
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class. The results are in perfect agreement with a perturbated conformal field theory scenario 
proposed by Cardy ^j. 
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A central element of the modern theory of bulk critical phenomena is the division into 
(bulk) universality classes. As is well known, the critical behavior near boundaries normally 
differs from the bulk behavior. In general, each bulk universality class of critical phenomena 
splits into several surface universality classes. 

The criteria by which the critical systems can be classified into different universality 
classes is a problem of much academic interest. Two-dimensional critical systems are pa- 
rameterized by the conformal anomaly c which is the central charge in the Virasoro algebra 
[sl, 14]. The conformal anomaly c can be obtained directly from the finite size corrections 
to the free energy for a two-dimensional classical system on infinitely long but finitely wide 
strip at a conformally invariant critical point. 

The asymptotic finite-size scaling behavior of the critical free energy fj^ and the crit- 
ical inverse correlation length C,i associated with the two-point correlations of in a two- 
dimensional spin system on a strip with infinite length and a width of Ai lattice spacings 
has the form 

lim M\fM - /oo) - 2Mfsurf = A, (1) 

hm A^C' = A, (2) 

M—^co 

where /oo is the bulk free energy, fsurf is the surface free energy and A and Di are the 
universal constants, but may depend on the boundary conditions (BCs). The index i distin- 
guishes different correlation lengths: for example, i = s for the spin-spin correlation length 
or i = e for the energy-energy correlation length. In some two-dimensional geometries, the 



values of A and Dj is known 
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6|, to be related to the conformal anomaly number (c), the 
highest conformal weight (A), and the scaling dimensions of the i-th scaling field (xi) of the 
theory 

Di = 2TTXi for periodic or antiperiodic BCs, (3) 

Di = TTXi for free, fixed and mixed BCs, (4) 

The principle of unitarity of the underlying field theory restricts through the Kac formula 
the possible values of c and for each value of c only permits a finite number of possible 
values of A. For the 2D Ising model, we have c = 1/2 and the only possible values are 
A = 0, 1/16, 1/2. The highest conformal weight A, and the scaling dimension Xi depends 
on the BCs. For the Ising model on infinitely long cylinder there are two different boundary 
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universal classes: periodic and antiperiodic with 

A = 0, Xs = -, Xe = I for periodic BCs, (5) 
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A = — , Xs = -, Xe = I for antiperiodic BCs. (6) 

For the Ising model on infinitely long strip there are three different boundary universal 
classes: free, fixed (H — ) and mixed with 



A = 


0, 


Xs — ^, 


Xe = 2 


for free and fixed (++) BCs, 
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A = 


1 

2' 


Xs = 2, 


Xe=2 


for fixed (+-) BCs, 


(8) 


A = 


1 
16' 


Xs -1-5 


Xe = 2 


for mixed BCs. 


(9) 



For fixed ++ (or H — ) boundary conditions the spins are fixed to the same (or opposite) val- 
ues on two sides of the strip. The mixed boundary conditions corresponds to free boundary 
conditions on one side of the strip, and fixed boundary conditions on the other. 

In the terminology of surface critical phenomena these three boundary universal classes: 
free, mixed and fixed (H — ) correspond to "ordinary", "special" and "extraordinary" surface 
critical behavior, respectively. 

Quite recently, Izmailian and Hu |7|, |8| studied the finite size correction terms for the 
free energy per spin and the inverse correlation length of critical two-dimensional (2D) 
Ising models on A^ x oo lattice and one- dimensional quantum spin model with periodic, 
antiperiodic and free BCs. They obtain analytic expressions for the finite-size correction 
coefficients a^, hk and c^ in the expansions 



M Um - foo) = 2fsurf + Yl J^V (10) 

k=l 

oo , 



fc=l 



C' = E^- (12) 

fc=l 

and find that although the finite-size correction coefficients ak, bk and Ck are not universal, 
the amplitude ratios for the coefficients of these series are universal and given by 

rs{k) = ^ = -j^^^-; r.(l)=3, r.(2) = ^, ... (13) 

ttk 2,"^ — 1 7 

-°W - 5 - (2^^- -\)i>. ^ ■^•W^^^- ^=P)^-H^. .... (14) 



for periodic BCs, where -B„ is the n-th BernouUi number {B2 = 1/6, B4 = —1/30, . . . 



^=«^^^^^^^^t^^ .„.-£, u.^f, ... (.5, 

r,{k) = ± = --— ; re(l) = -12, r,{2) = 120, .... (16) 
for antiperiodic BCs, 

r.ik)^'± 4. .._. .„. 240 



a, (22^-1 -l)52fc' 

,,, Cfc 4fc(32^-i + l) 
re(fc) = — 



s(l) = 24, 


rs{2) = 


7 ' ■■ 


• (17) 


e(l) = 96, 


r,(2) = 


-960, . . 


• • (18) 



a, (22*^-1-1)^2^.' 
for free BCs. 

In this paper we present exact calculations for a set of universal amplitude ratios for the 
two-dimensional (2D) Ising models on A^ x 00 lattice with the special boundary conditions 
studied by Brascamp and Kunz (BK) [9]. They considered a lattice with 2A/' sites in the 
X direction and Ai sites in the y direction. The boundary conditions are periodic in the x 
direction; in the y direction, the spins are up (+1) along the upper border of the resulting 
cylinder and have the alternative values along the lower border of the resulting cylinder. It 
was shown 10| that the asymptotic finite-size scaling behavior of the critical free energy fj^ 



of the Ising model on infinitely long strip with Brascamp-Kunz boundary condition has the 
form 

hm M\fM - /oo) - 2Mfsurf = -^, (19) 

which is consistent with the conformal field theory prediction for the mixed boundary con- 
dition (see Eqs. ([T]), (jl]) and (jH])) although the mixed boundary condition and the BK 
boundary condition are different on one side of the long strip. 

We obtain analytic equations for a^, bk and Ck in the expansions given by Eqs. (TTOj) . ( ITTj) 
and flT2|) and find that universal amplitude ratios the two-dimensional (2D) Ising models on 
A^ X 00 lattice with mixed BCs are given by 

r,{k) = ^ = -i^; r,(l) = -24, r,(2) = 240, . . . (20) 

cu 2'^''^^k 
r,{k) = ^ = -^; re(l) = -48, r,{2) = 1920, .... (21) 

O'k -D2fc 

As far as we know, no previous RG arguments, analytic calculations, or numerical studies 
predict the existence of this whole set of universal amplitude ratios. 



Consider an Ising ferromagnet on an N' x A4 lattice. The Hamiltonian of the system is 

|3H=-JY,s^Sj, (22) 

where /3 = {kBT)~^, the Ising spins Sj = ±1 are located at the sites of the lattice and the 
summation goes over all nearest-neighbor pairs of the lattice. There are a few boundary 
conditions for which the Ising model has been solved exactly. Among them is the special 
boundary conditions studied by Brascamp and Kunz (BK) [S] . We consider a transfer matrix 
acting along the A^ direction. If Aq, Ai and A2 are the largest, the second-largest and the 
third-largest eigenvalues of the transfer matrix, in the limit A/" — t- 00 the free energy per 
spin, fM = liniAr^oo F/2M, 



.v^oo 2^{M + 1) 
and the inverse spin-spin correlation length, ,^~^, and the inverse energy-energy correlation 

length, ^~^, are 

^•^ = >lVl^'''^°' C' = ln(Ao/Ai) and C' = In (A0/A2). (23) 

where F is the total free energy. The three leading eigenvalues of the transfer matrix (Aq, 
Ai and A2) can be obtained from exact expression for the partition function of the Ising 



I9 



model on an 2Af x Ai rectangular lattice under under Brascamp-Kunz BCs [9|]: 

2M+1 / 

Ao = C7,exp<(^y^c.,^^^^^^)^ (24) 



{, 2M+1 
m=0 

^^ - ^'-P {-^-" (aTTi) 4l/^ iliXTTT)) } ■ ^''^ 



2M+1 



where 



C = (V2e^]^ ( ^ Y 

'' V / V4 cosh (2Aru;^(0)) cosh (2Arw^(7r/2))y 

and a lattice dispersion relation uj^{x) is implicitly given by 



UJ^{x) = arcsinh Y 2 sinh /i -|- sin^ x 
with /i = I In sinh 2 J. At the critical point /i = /ic = (Jc = | In (1 + \/2)) one then obtains 

uJo{x) = arcsinhsinx 



and 

1 

-\M f 1 \ 477 



Co 



v^j 



4 cosh (2A/'arcsinhl)/ 

Then the critical free energy Jm-, critical spin-spin correlation length E,s and critical energy- 
energy correlation length ^e of Eq. ( l23l) can be written as 

'" = 2(Am)'"^-2(A?TT)'"<^ + ^' + 2(ATTT) 5/°(2(^)-P^) 



Using the Euler-Maclaurin summation formula ll| the asymptotic expansion of the critical 



free energy fj^ can be written in the following form 

(M + DUm- U = 2/.„; - ± ,;^^^,, i^rj^, ) , (30) 



^_„(2«:)!(2A + 2) \2{M + i: 

1 / \ 3 1 / \ 5 

TT i/7r\ i/TT 



where 



2/s«r/ 24(_/^ + 1) 2880 \M + l) 48384 \M + l' ^ 



1 2r^ 

/oo = -ln2 + — (31) 

Z vr 

/sur-/ = -^ln(2 + 2V2) (32) 



and A2fc is the coefficients in the the Taylor expansion of the uo{x 

-^2p 2p+l \ _ 1 \ _ 2 

(2^^^ ' Ao-l,A2--- 



oo 



^o(^) = E 7^ ^''■''' Ao = 1, A2 = -|, A4 = 4, ... (33) 

p=0 



Using the Taylor expansion of the 070(3;) given by Eq. ( l33l) the asymptotic expansion of the 
critical spin-spin correlation length ^^ and critical energy-energy correlation length C,e can 

be written as 

^^ ' = IJ (2^ [kmTT)) ' ^^^^ 

TT l/^7r\ I/ttA 



(2A;)! V-M + 1 



^- ' - 5Z 7^ T7TT ' (^^) 



fc=0 



27r 2/7r\ l/vr^ 



>1 + 1 3\M + lJ 3\M + 1 



+ 7T ^^^ + 



Equations (EU]), ([MD, and ([3S]) imply that the ratios of the amphtudes of the {M + 1) -(2*^+1) 
correction terms in the spin-spin correlation length, the energy-energy correlation lengths, 
and the free energy expansion, i.e. bk/dk and Ck/ak, should not depend in detail on the 
dispersion relation {uq{x)) as given by Eqs. fl20|) and fl2Tl) . 

To check the applicability of these results, we study the anisotropic Ising model with 
coupling constant J and 7J along the horizontal and vertical directions, respectively, with 
< 7 < 00. At the critical point Jc, where Jc is defined by sinh2JcSinh27Jc = 1, we obtain 
that the dispersion relation for the anisotropic Ising model Uq^^° ^°^ (x) is now given by 

Wg ^'"' (a;) = arcsinh (sinx sinh 27 Jc) (36) 

The asymptotic expansion of the critical free energy /x, the critical spin-spin correlation 
length ^s and critical energy-energy correlation length ^g is given by the first line in the Eqs. 
(!30|) . (IMIl . and (1351) . where X2k is now the coefficients in the the Taylor expansion of the 

1 . . 1 J. • (anysotrop) / \ 

dispersion relation cJq [x): 



^{anysotrop) /-^ ^ '^ ^2p ^2p+l 



00 , 



where 



Ao = sinh 27 Jc, A2 = — sinh 27 Jc cosh^ 27 Jc, 

o 



It is easy to see that the ratios hk/ak and Ck/ak does not depend in detail on the dispersion 
relation {ujq^^°^{x)) and Eqs. fl20l) and fl2T]) holds for all anisotropy 7. 

In [7|, we have shown that the Ising model on the square, honeycomb and plan-triangular 
lattices, and the quantum spin model have universal amplitude ratios, i.e. we confirmed that 
such models are in the same universality class. It is reasonable to expect that the ratios of 
Eqs. ( I20I) and ( 12T1) are valid for the same set of models with mixed BCs. 

The leading terms of Eqs. fl30l) . flMl) and fl35l) are consistent with Eqs. ([I]) - dl]), «.e. ai, 
hi and Ci are universal. Equations ([1]) and ([2]) implies immediately that their ratio is also 
universal, namely ^^(1) = Di/A and re(l) = D2/A, which is consistent with Eqs. fl20l) and 
(|2T1) for the case k = 1 

r,(l) = ^ = -2A (38) 

re(l) = ^ = -48 (39) 



The finite-size corrections to Eos. Q and ([2]) can be calculated by the means of a 



perturbated conformal field theory 2; 
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In general, any lattice Hamiltonian will contain 



correction terms to the critical Hamiltonian Hr 



j-M/2 
„ J-M/2 



(f)p{v)dv, 



(40) 



where Qp is a non- universal constant and 0p(f) is a perturbative conformal field. Below 
we will consider the case with only one perturbative conformal field, say (piiv). Then the 
eigenvalues of H are 

rM/2 

En = En,c + gi < n\(t)i{v)\n > dv + . . . , (41) 

J-M/2 

where En^c are the critical eigenvalues of H. The matrix element < n\(j)i{v)\n > can be com- 
muted in terms of the universal structure constants {Cnin) of the operator product expansion 
2|: < n\(j)i{v)\n >= (27t/M.)^^ Cnin-, where xi is the scaling dimension of the conformal field 
bi{v). The energy gaps {En — Eq) and the ground-state energy {Eq) can be written as 

Xl-l 



27r / 27r 

En- Eq = — Xn + 2ngi{Cnin - Coio) ( -r- 

(r) \ Xl-l 



+ ... 



(42) 
(43) 



Note, that the ground state energy Eq, the first energy gap {Ei — Eq) and the second energy 
gap {E2 — Eq) of a quantum spin chain are, respectively, the quantum analogies of the 
free energy f{N), inverse spin-spin correlation length ^~^(A^), and inverse energy-energy 
correlation length ^~^(A^) for the Ising model; that is, 

MfM^-Eo, C\N) ^ El - Eo = As, and C\N) ^ E2 - Eo = A,. (44) 



For the 2D Ising model, one finds 13| that the leading finite-size corrections (1/A^'^) can be 
described by the Hamiltonian given by Eq. (140|) with a single perturbative conformal field 
(f)i{v) = L^2(^) with scaling dimension xi = A . 

In order to obtain the corrections we need the matrix elements < n|L^2(^)l^ >5 which 



have already been computed by Reinicke 14j ]. 



<A + r\L^_,\A + r> = (^] 



11520^ ^V 24 (A + 1)(2A + 1) 



(45) 



The universal structure constants C212, Cm and Com can be obtained from the matrix element 
< n\L'^2{'^)\^ >= {27T/AiY' Cnini where x/ = 4 is the scaling dimension of the conformal 
field Ll^^v). 

At the critical point the spectra of the Hamiltonian with free, fixed and mixed BCs are 
built by the irreducible representation A of a single Virasoro algebras with possible values 
of A are 0, |, ^. We denote by A the highest weight, and by A + r, the r-th level having 
degeneracy c/(A, r) of irreducible representation of the Virasoro algebra. A state will be 
labeled by \n >~ | A + r >. 

For mixed BCs the ground state 10 >, first excited state 11 >, and second excited state 
1 2 > are given by [2, Il5l |: 



|0> = |A = ^,r = 0>, (46) 

|1> = |A = l,r = l>, (47) 

|2> = |A = 4,^ = 2 >. (48) 

lb 

After reaching this point, one can easily compute the universal structure constants 6*2^25 
Cm and Cqiq for mixed boundary conditions. The values of Cqiq, Cm, C212 can be obtained 
from Eqs. ( I45ll - ( l48l) and given by: 

Cqio = -7/1440, Cm = 1673/1440, C212 = 13433/1440. (49) 

Equations (H2l) and ( l43l) implies that the ratios of first-order corrections amplitudes for 
En — Eq {in^) and —Eq {j'm) is universal and equal to (Cow — C'nin)/C'o;o; which is consistent 
with Eq. fl20l) for the case n = l,k = 2 

r,(2) = ^°'° ~ ^^'^ = 240 for mixed BCs (50) 

Cqio 

and with Eq. ( l2Tl) for the case n = 2,k = 2 

r,(2) = '^o'Q ~ ^^^^ = 1920 for mixed BCs (51) 

Cqio 

In this paper we present exact calculations for a set of universal amplitude ratios for the 
two-dimensional (2D) Ising models on A^ x 00 lattice for mixed BCs universality class. We 
find that such result are in perfect agreement with a perturbated conformal field theory 



scenario proposed by Cardy 
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